Decoherence by electromagnetic fluctuations in double-quantum-dot charge qubits 
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We discuss decoherence due to electromagnetic fluctuations in charge qubits formed by two lateral quantum 
dots. We use effective circuit model to evaluate correlations of voltage fluctuations in the qubit setup. These 
correlations allows us to estimate energy (Ti) and phase (T2) relaxation times of the qubit system. Our theoret- 
ical estimate of the quality factor due to dephasing by electromagnetic fluctuations yields values much higher 
than those found in recent experiments, indicating that other sources of decoherence play a dominant role. 
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I. INTRODUCTION 

Solid-state semiconductor lateral quantum dots are strong 
candidates for the physical realization of qubits. These artifi- 
cial systems can be designed to allow for the observation of 
coherent oscillations between their quantum states. Since its 
first proposals,' ^ a wide variety of experiments have demon- 
strated control over the spin degree of freedom of confined 
electrons in quantum dots,"* "* as well as charge states.^ ** Solid- 
state quantum computer architectures with qubits encoded 
in dopant atoms in semiconductor crystals have also been 
proposed."^ Quantum dots present the ubiquitous advantages 
of being manufactured from highly developed semiconduc- 
tor technology and may offer easier scalability, the latter be- 
ing key in enabling the manufacturing of large-scale quantum 
computers in the future. A drawback to their use in quantum 
computers is that they also couple rather effectively to external 
degrees of freedom which lead to decoherence. 

Semiconductor qubits are susceptible to various decoher- 
ence mechanisms. Hyperfine coupling to lattice n uclea r spins 
reduces the phase coherence of electron spins, '°" while 
qubits based on the charge degree of freedom are particu- 
larly sensitive to decoherence mechanisms related to charge 
motion, such as coupling to phonon modes and to charge 
traps in the substrate. Several of these sources have been 
investigated-'i^"^ 

In this paper we focus on charge-based qubits. The simplest 
realization of a charge qubit is a doubl-quantum-dot (DQD) 
system with an odd number of electrons, as shown schemat- 
ically in Fig. [T] So far measurements of quality (Q) factors 
of coherent oscillations in these systems have yielded rather 
low values in the range of 3-9.* ^ In an effort to identify the 
main sources of decoherence, theoretical estimates of the Q 
factor have been carried out assuming mainly the coupling to 
acoustic phonons.t^^lEilHowever, a discrepancy of at least one 
order of magnitude remains between the experimental value 
and the theoretical estimates, with the latter predicting larger 
Q factors. This discrepancy indicates that the phonons may 
not be the dominant noise source in current experimental se- 
tups. Thus, an investigation of other possible environmental 
decoherence mechanisms is in order. Here, we consider the 
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FIG. 1 . (Color online) Schematic representation of a double quantum 
dot setup. 



coupling of the DQD charge-based qubit systems to voltage 
fluctuations in the gates. 

This paper is organized as follows. In Sec. Ill] we define the 
DQD effective Hamiltonian and the interaction between the 
qubit and the gate-voltage fluctuations. In Sec. [Ill] we intro- 
duce the effective circuit model that describes the DQD and 
the electromagnetic environment, as well as the Hamiltonian 
of the latter. In Sec 
rameters and in Sec 



Vlwe estimate the equivalent circuit pa- 
VI] we calculate upper bounds to the de- 
coherence rates and Q factors. Our main finding is that volt- 
age fluctuations cause only a very small decoherence effect 
in DQD charge qubits. Since double-dot spin qubits are also 
susceptible to decoherence due to charge motion, in Sec. |VII| 
we use results from our circuit model of electromagnetic fluc- 
tuations to estimate decoherence rates for those system. Our 
conclusions are presented in Sec. |VIIl] 



II. HAMILTONIAN OF THE DOUBLE-QUANTUM-DOT 
SYSTEM 



The Hamiltonian of a DQD can be separated into a quantum 
part related to the occupation of energy levels on each dot and 



a classical part that quantifies the charging energy, 

H = ^einc\^Cin+^e2ncl„C2n+E{Ni,N2), (1) 
n n 

where c]^„ and Ci„ are creation and annihilation operators of 
the state with energy Sin in the left {i — 1) or right dot {i — 2). 
The dot occupation numbers are defined as TV^ = ^^ c^in^in 
while the total charging energy is given b}P2 
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The capacitances and voltages shown in Eqs. (|2|-(|5]l are de- 
fined in Fig. I2] C1.2 is the sum of all capacitances attached to 
dot 1 or 2: ct,2 = Cria + C3i,2 + C„. 

For the purpose of our analysis, the Hamiltonian can be 
greatly simplified. Notice that the DQD qubit can be viewed 
as a double-well potential where an unpaired electron oscil- 
lates between both quantum dots by tunneling through the po- 
tential barrier. Spin degrees of freedom can be neglected. By 
adjusting the gate voltages, one can set the system near the 
degeneracy point E{1, 0) = £'(0, 1), in which case the logi- 
cal states of the qubit correspond to the electron being on the 
left or right, \L) {Ni = 1 and N2 = 0) and \R) (Ni = 
and N2 = 1), respectively. The typical single-particle level 
spacing within each quantum dot is assumed sufficiently large 
so that only one level on each dot needs to be considered at 
low enough temperatures. The barrier height A determines 
the tunneling rate between the dots and can be adjusted by a 
gate voltage while a bias e between the two dots can also be 
applied through two independent plunger gate voltages. The 
dynamics in the DQD qubit is then governed by the reduced 
two-level Hamiltonian 

Hs='- {\L){L\ \R){R\) + I {\L){R\ + \R){L\) , (6) 

with the constraint that \L){L\ + \R){R\ = 1. The fields e 
and A represent the interdot bias and the interdot capacitive 
coupling, respectively. 



Electromagnetic noise is introduced into the DQD qubit 
system by means of gate voltage fluctuations. These fluctu- 
ations may originate from the voltage sources and the ther- 
mal noise in the transmission lines, and introduce decoherence 
into the qubit system through interactions with the electrons 
in the quantum dots. While the former can be substantially 
reduced by careful filtering, the latter is less controlled. Here 
we will focus on the noise coming from the plunger gates. The 
effect of voltage fluctuations in the gate electrodes is captured 
by the qubit-environment interaction 



HsB = e7]iSVgi - SVg2) m{L\ - \R){R\) 
where rj is the capacitive lever arm coefficient. 
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Depending on the particular qubit setup, other sources of elec- 
tromagnetic noise may also exist, such as bias and current- 
voltage fluctuations. They can affect not only the qubit coher- 
ent dynamics but also the state measurement. For the sake 
of maintaining some generality in our study, we will only 
treat electromagnetic fluctuations which can be expressed as 
Eq. dTli. In addition, we will model the voltage fluctua- 
tions through frequency-dependent impedances along the gate 
transmission lines. 



III. HAMILTONIAN FOR THE ELECTROMAGNETIC 
ENVIRONMENT 

The effective circuit of a double quantum dot setup is shown 
in Fig. |2] The effect of the electromagnetic environment is 
modeled by the frequency-dependent impedances Zi 2(w). In 
the experimental setups, the voltage lines typically run parallel 
to each other over several microns or more. In order to take 
into account any capacitive coupling between the lines, we 
introduced capacitance C12 into the circuit. 
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FIG. 2. Circuit representation of a double quantum dot system cou- 
pled to an electromagnetic environment through metallic gate elec- 
trodes. Source and drain electrodes are assumed grounded. The 
quantum dots are denoted by QDi and QD2. 



The impedances Zi,2(w) can be modeled by means of a 
transmission line with distributed elements, which stems from 
the fact that the source of noise in our circuit is spatially 
distributed along a finite length. Let us consider first each 
transmission line independently, as shown in Fig. [3] whose 
impedance Zi (w) can be represented by an infinite ladder net- 
work of identical inductors Ln and capacitors Cu (see Ref. 
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Typically, it would be necessary to estimate the values of the 
spatially distributed resistance, capacitance, and inductance 
in the circuit, but the choice to model the impedance as a 
LC transmission line can be made because it is known that 
through a (not necessarily trivial) normal-mode transforma- 
tion, any RLC or RC transmission line can be written as an 
infinite LC ladder network. The elements Ch and Lu of the 
transmission line can be determined from two real parame- 
ters of the real: the cutoff frequency lOc and the low frequency 
asymptotic limit to the characteristic impedance Z{uj = 0). In 
an semi-infinite line, Re{Z(a;)} = when uj > lOc- Hence, 



(10) 
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Zi{uj = 0), on the other hand, can be calculated by taking the 
low frequency asymptotic limit of Eq. (|9]l. It is straightfor- 
ward to see that this limit yields 



Adding the capacitive coupling between the voltage trans- 
mission lines, we obtain the following environmental noise 
Hamiltonian: 
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The cross term complicates the task of finding the normal 
models of the environment and an alternative approach was 
adopted. 



IV. DOUBLE-DOT JUNCTION 

The double junction solution is based on the original solu- 
tion for a single-dot junction treated in detail by Ingold and 
Nazarov.™ The nontrivial aspect of our extension of the cal- 
culations in Ref. "24" is the inclusion of the gate capacitances 
(see Fig. pi. 

We start with the setup shown in Fig. |4] Following a 
straightforward application of Kirchhoff's laws, we find the 
relations 



with 



Vl = (il + il2 + irn) Zi + Vgi, 
V2 = {il - il2 - im) Z2 + Vg2, 



Vgl - Vg2 = il2 Zi2 



(14) 
(15) 



(16) 



(11) 



where R is an ohmic resistance. 

To introduce noise, the transmission line is decomposed 
into normal modes. In quantized form, the charge at the Ith 
node, Qii, and the flux (/); i are conjugated variables obeying 
the commutation relation [</>;, i, Qr.i'] — ieSi^i'Si^r. Following 
the standard procedure, we define the Hamiltonian governing 
the flux and charge fluctuations along such transmission as 



Htj = 
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Notice that Cgi represents the capacitive coupling between the 
quantum dots and their respective gates, while Cu and L^ 
represent the capacitive and inductive terms, respectively, at 
each rung in the transmission line. 
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and 

Ui — U2 = im Zm , 

where Z12 — (ia;Ci2) and Z„i — {iujCm)~ ■ 



(17) 




FIG. 4. Circuit of a double-dot junction system coupled to two volt- 
age sources through noisy lines. 

We begin by eliminating ii2 and i„i in Eqs. ([14]) and ( [T5| ) 
with the help of Eqs. (16 1 and (17 1, and proceed to write Vgi 



and Vg2 in terms of Ui and 1/2- For this purpose, we notice 
that 



FIG. 3. Circuit representation of the electromagnetic environment as 
a transmission line. 



Vgl = C/i + {il +i,n) Zgi, 

Vg2 ^ U2 + («2 - im) Zg2, 



Ci Ui = Qi, and ii — Qi, with i = 1, 2. Eliminating Vgi, Vg2 
and rewriting Vi^ V2 in terms of ii,i2, Qi and (52, we obtain, 
in matrix notation 



where 
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with 
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The matrix C is defined as 
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Notice that when we set Cm — and C12 = in Eq. ( [T8| , 
we decouple the two halves of the circuit and obtain two inde- 
pendent equations for each half of the circuit. 

In analogy to the single dot-junction circuit, the Hamilto- 
nian for the environment in this case can be written as 
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We find that 
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where a;„j = l/\/I~C^, and 
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(27) 



C„i(i) — — 



ei^ 



sm{ujnit) V'miO) 



+ COs{uj„it) ipni{0) 



(28) 



In addition, it is easy to show that the relation Qi 



~ X]^i irii also hold. Thus, we can write 

Q, (t) + - / dt' y, (t - i') <^g, (i') = /^, (t) (29) 

e Jo 

where the parameters {C„i, L„i} must be chosen such that 

COs{uJntt) .., . 1 
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A. Voltage correlation functions 

We now turn the charge equation of motion, Eq. (J29]l, into 
one for phase fluctuations, by using the following relationship: 



"Pg 



{Vg - V) 



(32) 



and 
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Substituting Eq. ( 33 1 into Eq. ( 32 1, we arrive at 



^g 



eQ 

he' 



(33) 



(34) 



which allows us to retrieve the phase fluctuation equation-of- 
motion. 

~ ^ /■* ^ e 

C-^g{t)+J dt'y{t-t')^g{t')^-lN{t). (35) 

Since we are interested in the behavior of a double-dot junc- 
tion, with each dot possessing its own charge and phase fluc- 
tuations, we will from now on represent these quantities in a 



vector notation, as seen in Eq. ( 35 1. By applying a Fourier 



transformation and substituting the random internal currents 
by external ones, we get 



iuZ^ ^[u) ■ (pg{uj) = jTpcrt (w), 



where 



Zt\u:) 



icoC. 



(36) 



(37) 



Now, substituting the external currents by appropriate gener- 
alized force matrix. 



we obtain 



Tg{uj) = - Xpert (w), 



iPg{uj) ^X^^{UJ) • Jg(w), 
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e\2 RejZtH} 



(38) 



(39) 



where the dynamical susceptibility matrix is given by 

KM-{::i)iu^\ (40) 

whose imaginary part is given by 



(41) 



Assuming that both transmission lines are at the same tem- 
perature, the generalized form of the fluctuation-dissipation 
theorem reads 






Hence, 
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We now turn to the fluctuations of the voltage at the dots. 
Since 



S)^-'-(t)^''^-'^'^-(l::>'<«' 



where 



we find the matrix equation 



Ci 

c. 



U = in^M, 

\ _ g — phLO ' 



where 



(45) 



(46) 



M^C-^ -C- Re {Zt{uj)} ■ & ■ {C-^)^ , (47) 



and 
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V. ESTIMATE OF CIRCUIT PARAMETERS 

We now proceed to make realistic estimates of the effec- 
tive circuit parameters. The double-dot system is maintained 
at very low temperatures, in the tens of millikelvin.^^ Typi- 
cally, ksT <C AE, Eci, Ec2, where AE is the mean level 
spacing in the dots. The wires leading to the double quantum 
dot are thermally anchored to a fridge at several temperature 
stages (4 K, 1 K, 100 mK, and 10 mK). The transmission line 
resistance Rj^ is estimated to be 50 fl for low temperatures (at 
or below 4 K) inside the dilution refrigerator, or 250 fl in the 



copper leads residing at room temperature P^ 

The resistance of the two-dimensional electron gas (2DEG) 
can be calculated using Drude's theory. ^^ The typical electron 
density in a high-mobility GaAs 2DEG is approximately n = 
10^^ cm^^, which leads to an average Fermi velocity of about 
vf — 10^ m/s. At subkelvin temperatures, mean-free paths 
in the 2DEG range from a few to up to one hundred fim?^ 
Choosing I — 10 /im, we arrive at a relaxation time r — 
l/vp ~ 100 ps, leading to an estimate of the low-temperature 
conductivity of 



ne^T 



4.2 X 10" 



S, 



with m* = 0.067me = 0.61 x lO^'^^ kg being the electron 
effective mass in GaAs. To calculate the resistance, we con- 
sidered a length I — 10 fim and a width w — 2.5 /im, yielding 
a sheet resistance for the 2DEG underneath the gate electrodes 

i? = p- ~ 95 r^/D, 
w 

where p = 1/cr is the resistivity of the 2DEG. This resistance 
is responsible for a dissipative drag effect,^^ that, for the sake 
of simplicity, will not be considered in our model. 

There is still one resistance left to be determined, which is 
the resistance of the metallic electrodes. This resistance can 
be determined by 



R = p 



I 
he' 



(49) 



where p is the resistivity of the electrodes, approximately 
0.022 X 10^^ Viva for a Au electrode at low temperature (< 
4 K). If we consider the electrodes to have a 10 /im length and 
a 30 run x 60 nm cross section, we can estimate the electrode 
resistance to be around 1 il, a small value that will also not be 
considered in our model. 

The capacitance C between the transmission line and the 
2DEG was estimated by solving the electromagnetic problem 
of a cylindrical conducting wire of radius r = 20 nm placed at 
a distance of d = 100 nm from an infinite grounded conduct- 
ing plate. Using the method of images, we can estimate the 
total electric potential of this system by integrating the elec- 
tric field along the line connecting the centers of the real and 
the image wires. This results in a capacitance per unit length 
of 25 aFZ/xm, and a total capacitance of 250 aF for a wire of 
10 pm. in length. 

Any inductive couplings along our voltage lines can be esti- 
mated as follows. For a metal electrode with rectangular cross 
section, the self-inductance in H/m is approximated a^ 



-t^rod '^ 2t 
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(50) 



where e is the aspect ratio of the electrode. For an electrode 
with an aspect ratio of 2, this equation yields i w 1 pH//im. 
Thus, a 10 pm long electrode gives us an inductance of 10 
pH. The parameters C = 250 aF and L — 10 pH, though 
useful as rough estimates to characterize circuits, will not be 



used in our model since they are very specific to the given cir- 
cuit. In fact, in order to estimate these circuit elements more 
precisely, more physical parameters of the circuit in question 
would be necessary. To determine the transmission line pa- 
rameters in our model, we will make use of Eqs. ([T0]l and 
( [TT| ) from Sec. |lIl]to give us a more general approach where 
we can model any transmission line given these two operating 
parameters. To give us a large enough window to operate our 
qubits, we set our cutoff frequency to Wc = 200 x 10^ rad/s. 
Table|l]summarizes the transmission line parameters that fully 
describe Zi{uj). 



TABLE I. Estimates for the transmission line parameters. 



Transmission line parameters 


Length I 


10 /im 


Transmission line capacitance Ct 


10 pF 


Transmission line inductance Lt 


10 pH 


Cutoff frequency lOc 


200 x 10^ rad/s 


Z{u = Q) = R 


1 Q, 



The gate capacitance Cgi 
is given by 



1, 2) for each quantum dot 



^'' ~ AK. 



(51) 



If we consider AK, « 4.5 mVpHUwe find C„, « 40 aF 



' g% 



9« 



Finally, we now estimate the tunneling parameters between 
the quantum dots and the 2DEG. These are given by a tun- 
neling junction with an impedance Zt = Rt + jXcj.- We 
can obtain a lower bound for the tunneling resistance Rt by 
estimating the inverse of the Coulomb blockade peak conduc- 
tance. In the regime F <C ksT, Gmax is given b>3il 



Gr, 



FT'' 



ikeT r' + F'- 



(52) 



where tunneling rates of an electron through the potential 
barrier into (or out of) each dot are assumed equal for the 
sake of simplicity (F' = F'') For an electron temperature 
in the dot T « 150 mK and a peak conductance height of 
2 X lO^'^e^ /hp^we find the tunneling resistance to be larger 
than or on the order of 10 Mfl. We can estimate the tunneling 
capacitance indirectly. We know the expression for the total 
capacitance of a flat disk to be 



C^ — 8€j-cqR. 



(53) 



Assuming i? ~ 80 nm as the radius of the quantum dot and 
£,. w 11 for GaAs at high frequencies, yielding a total capaci- 
tance Ci w 60 aF for each quantum dot. 

From the total capacitance we can estimate the interdot ca- 
pacitance between dots 1 and 2, since 



G„ 






(54) 



where i ^ j. For AV™ w 0.4 mVp^'^we find G™ « 5 aF. 



The total capacitance for each quantum dot, as seen previ- 
ously, is the sum of all capacitances attached to the dot. As 
such, by knowing C„j = 5 aF and Cxi ~ 40 aF, we find Cgi 
« 15 aF 

Using these estimates for the circuit elements, we are able 
to determine the distributed parameters of our noisy transmis- 
sion lines. According to Eq. ( [TT] i, if we assume a cutoff fre- 
quency of Wc ^ 10^^ Hz, we find Lu ^ 1 pH//xm and Cu ^ 1 
pF//im. 

In principle, one can also consider the ground (2DEG) to 
be a source of noise, and as such it can also be modeled by 
means of a frequency-dependent impedance. This would re- 
quire however an appropriate estimate of the inductance along 
the 2DEG. We did not carry out such an estimate. How- 
ever, we attempt to take into account the coupling between 
the quantum dot leads. This coupling is given by the lumped 
capacitance C12, as shown in Fig. [2] This capacitance was 
estimated to be approximately 20 aF by means of numerical 
multipole expansion calculations performed by a field solver 
software.'22l We summarize in Table HH the relevant circuit pa- 
rameters necessary to fully characterize the DQD setup. 



TABLE II. Estimates for the DQD circuit parameters, i — 1,2, cor- 
responding to each quantum dot 



Circuit parameters 


Transmission line capacitance Cu 


1 pF//im 


Transmission line inductance La 


1 pHJfini 


Interdot capacitance Cm 


5aF 


Tunneling capacitance Cn 


40 aF 


Tunneling resistance Rn 


> 10 Mt^ 


Gate capacitance Cgi 


15 aF 


Total quantum dot capacitance d 


60 aF 


Capacitive coupling between transmission lines C12 


~20aF 


Electrode resistance Ri 


1 n 



VI. BOUNDS ON DECOHERENCE RATES AND Q 
FACTORS 

Through the fluctuation-dissipation theorem, we can re- 
late the impedance Z\^i{lS) to a source of electromagnetic 
gate fluctuations 5V^i.2. These gate fluctuations Wgi = 
Qo,i/Cgi{i ~ 1, 2) can be determined through the diagonal- 
ization of the Hamiltonian in Eq. (fTSll. We consider in this 
paper the case of Johnson-Nyquist noise. Following the 
standard procedure, we relate the energy relaxation rate to the 
power spectrum of voltage fluctuations, 



71 



sin ry 
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where 



Ssui'^) 



Ssu{^/h) + Ssui-^/n) 



\[6Ui{lu) - 6U2{l^) 



(55) 



(56) 



and tan 77 = A/e, see Eq. (46 1. Using Eq. ^, we obtain 
sin^ 77 A 
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■ coth 



A 
2kBT 



M 



(57) 



where 



M = Mil + M22 - M12 - M2i- 



(58) 



and Rk is the resistance quantum (= h/e^ ~ 25.8kri). From 
this expression we can calculate the energy relaxation and de- 
phasing times. 



Ti = 1/71 



HRk tanh(A/2fcBT) 



sin rj A 



and 



T2 = 
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M 
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Hereafter, we will assume zero bias (e — 0), in which case 
rj = 7r/2 and T2 = 2Ti. The quality factor of the quantum 
oscillations is then given by 



Q — ^osc T2 — 



7r7i 



(61) 



where Wosc is the frequency of quantum oscillations observed 
in the DQD system, as defined bjCS 



Wo 




(62) 



with A being the potential barrier height between quantum 
dots, as shown in Eq. (|6]), and 72 being defined as 



72 



^ iyi2Ay) coth ' ^ 



(63) 



where 1/ is the bath spectral function, defined as 

H^) = ~{{\SUiiu;)f) + {\SU2iu;)f 

- {dU^{uj)6U2{io)) - {6U;{u)SUi{uj)) }.(64) 

The operating frequency uj — 2A/h is fed to the circuit by 
the voltage generators and carried through the gates. The other 
terms in Eq. ( |62| ), as it turns out, are small enough corrections 
to the operating frequency so that they may be ignored. Thus, 
from now on we will assume cJqsc = ^■ 

We will now analyze in detail two different scenarios: one 
where the transmission lines are decoupled, while the other 
includes the capacitive coupling C12 between transmission 
lines, as seen in Fig. l2] 

A. Case (i): Decoupled transmission lines 

It is useful to look at the case where there is no coupling 
between the electrodes. The decoherence introduced by the 
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electromagnetic voltage fluctuations can still be analyzed us- 
ing Eqs. (|57]i-(|6T]l, but some simplifications to the impedance 
matrix are now possible. This case corresponds to having 
Ci2 = 0, so the matrix Z from Eq. ([T9]l is reduced to 



zA\ 






-z. 



Qrn 
C2 

C2 



(65) 



In this case we observe the highest possible quality factors for 
our double-quantum-dot setup, as seen in Figs. [5] and l6] 



3. 
O 




^/[Hz] 



FIG. 5. Qubit quality factor as a function of frequency for two de- 
coupled semi-infinite transmission lines, with temperature T — 150 
mK and the circuit parameters presented in Table [HI 



where Ri — Z{uj = 0) = \/Lti/Cu, as reported earlier, and 
with the assumption that Ci = C2, Cgi = Cg2, and Ri = i?2- 
This, combined with the fact that Re {Zt{u})} — >■ Re {Z{uj)} 
for w — > 0, yields 



lim Q{io) = 



.9 X 10 



-7r 



']) 



V, 



(67) 



where we notice a linear dependence of Q with respect to v, 
as can also be evidenced in the log-log graph shown in Fig. 
[7] While Ri is an important modeling parameter for the trans- 
mission lines, it is also clear that Cu and Lu ultimately influ- 
ence how quickly this linear regime establishes itself once we 
move to lower frequencies. 




[Hz] 



FIG. 7. Quality factor as a function of frequency for two decoupled 
transmission lines represented in a logarithmic scale with the same 
parameter values as in Fig. [5] 




[Hz] 



FIG. 6. Qubit quality factor as a function of frequency for v < 20 
GHz and two decoupled semi-infinite transmission lines with the 
same parameter values as in Fig. B] 

If we look back at Eq. ( |6T| and take its asymptotic limit for 
low frequencies, Eq. (|65|l is then reduced to 



Ri[l 



Ci 



-i?l 



Cm 
C2 



^■^2 7^ 



i?: 



('-«?);• 



(66) 



Turning our attention now to higher frequencies, we notice 
an important characteristic of the transmission lines. The real 
part of the transmission line impedance Re {Z(w)} has a cut- 
off frequency given by v^ — uJc/'^''^- In Fig-[8l it can be seen 
that as w -> lOc, Re {Z(a;)} — > 0, making Re {Zt(w)} — > 
as well, causing the quality factor Q to diverge at w = oJc- 



B. Case (ii): Capacitively coupled transmission lines 

Inserting now the inter-capacitive coupling C12 estimated 
in Sec. W\ we obtain the quality factor Q as a function of 
frequency v shown in Figs. |9j [TO] and [TT] In Fig. |9] we 
can clearly observe the quality factor diverge at the frequency 
Vc — 320 GHz, corresponding to the cutoff frequency. From 
now on we shall restrict our discussion to operating frequen- 
cies under 20 GHz (Fig.fTT), which are more realistic for prac- 
tical implementations of qubit operations. 

It is interesting to observe the influence of temperature on 
the decoherence introduced into the system by voltage fluc- 
tuations. We show below, in Figs. 12 and 13 a family of Q 
factor curves as a function of operating frequency v for tem- 
peratures ranging from 50 mK all the way to room temper- 
ature. As temperature increases, more environmental modes 




FIG. 8. Real part of the impedance Z(u}) as a function of tfie fre- 
quency u. Transmission line parameters are defined in Table [ll| 




[Hz] 



FIG. 10. Qubit quality factor as a function of frequency represented 
in a logarithmic scale. The circuit parameters utilized are the same 
as in Fig. [9] 



3 




[Hz] 



FIG. 9. Qubit quality factor as a function of frequency, with temper- 
ature T = 150 mK. The circuit parameters utilized are presented in 
Table In] 



are available for the system to couple with, effectively increas- 
ing dissipative effects. 

We can also observe the influence of the inter-capacitive 
coupling Ci2 on the quality factor, as seen in Figs. 14 and 
15 For weaker coupling, i.e., smaller C12, the quality factors 



are higher, as C12 approaches the limiting case of decoupled 
lines. Note that Q{lo) will still not reach the same levels of the 
decoupled case due to the presence of the capacitance Cm- 

In Tables [III] and |IV| we present the results of calculations 
for the decoherence time T2 and the Q factor for several dif- 
ferent values of temperature T and inter-capacitive coupling 
C12. It is easy to understand why higher temperatures degrade 
decoherence times in qubit operations. We can consider two 
extreme cases, namely, one where the electrical leads are in- 
side a dilution refrigerator and another where they are at room 
temperature. We will also consider an operating frequency 
V = Lo /2tt of 10 GHz. First, let us assume that leads con- 
nected to the gate electrode are inside the dilution refrigerator. 
In this case, a temperature T — 150 mK results in a relaxation 




i^[Hz] 



FIG. 11. Quality factor as a function of frequency for i^ < 20 GHz 
and two decoupled semi-infinite transmission lines with the same pa- 
rameter values as in Fig. [9] 



time Ti — 88 ns and a decoherence time of T2 — 176 ns. This 
scenario yields a quality factor of Q w 1760. If we consider 
now the case where the leads are at room temperature, we esti- 
mate the relaxation time and the dephasing time to be approx- 
imately 76 ps and 152 ps, respectively, resulting in a quality 
factor of Q — 1.5, more than 1000 times lower. A much more 
interesting analysis stems from varying the inter-capacitive 
coupling between the transmission lines. For higher values 
of C12, it would be intuitive to expect both transmission lines 
to be more strongly coupled, meaning that decoherence in 
the system would be weaker since voltage fluctuations in the 
two lines would be correlated. As it turns out, however, the 
stronger coupling between transmission lines results in larger 
off-diagonal terms in the matrix of voltage correlations de- 
fined in Eq. ( |46| l. If we look at Eq. ( |6T] i once more, it is 
easy to see that larger off-diagonal terms subtracted from the 
main diagonal correlation terms results in smaller Q factors, 
as evidenced by the behavior of the family of Q factor curves 
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FIG. 12. Qubit quality factor as a function of operating frequency 
for temperatures T = 50, 150, 250, 500 mK, and 1 K. The circuit 
parameters utilized are presented in Table [HI 



FIG. 14. Quality factor as a function of operating frequency for tem- 
perature T = 150 mK and inter-capacitive couplings C12 =0, 1.3 
, 10, 20, and 50 aF The circuit parameters utilized are presented in 
Table ini 
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FIG. 13. Logarithmic representation of the qubit quality fac- 
tor as a function of operating frequency for temperatures T = 
50, 150, 250, 500 mK, and 1 K. The circuit parameters utilized are 
presented in Table [ll| 



3i 
o 











C12 

— 
-- 1 aF 

— 10 aF 
20 aF 

— 5CaF 








^ 


^ 


^: 


^ 


^ 







[Hz] 



FIG. 15. Logarithmic representation of the qubit factor as a func- 
tion of operating frequency for temperature T = 150 mK and inter- 
capacitive couplings C12 = 0, 1.3, 10, 20, and 50 aF The circuit 
parameters utilized are the same as in Fig. [14] 



in Fig. [T4|for different values of inter-capacitive coupling and 



the calculated values presented in Table IV 



VII. 



ELECTROMAGNETIC FLUCTUATIONS IN 
DOUBLE-DOT SPIN QUBITS 



Decoherence due to the coupling between orbital (charge) 
and environmental degrees of freedom also occurs in certain 
spin-based quantum dot qubits. For instance, in the double- 
dot system introduced by the Harvard group^, the computa- 
tional basis is formed by the singlet and the 5^ = triplet 
states of a DQD system possessing an overall excess of two 
electrons. Single qubit operations are performed by modu- 
lating the gate-voltage difference between the two dots, as 
well as through the coupling to an inhomogeneous Overhauser 
field. For instance, calling the singlet states "up" and triplet 



state "down" pseudospin states, we can write the following 
pseudospin Hamiltonian 



iis = 



dJ 

Ik 



(e) 



67? {Vi-V2),(T,+HhF, 



(68) 



where Hhf describes the coupling to the Overhauser field 
and J(e) is the effective exchange coupling. The latter can be 
calculated in second-order perturbation theory, 



J(e) = At^ 



1 



1 



(U-S) ([/ + £) 



(69) 



with t = A/2 denoting the interdot tunneling matrix element 
and U = Eci = Ec2 representing the dot charging energy 
(for the sake of simplicity, we assume equal charging energy 
for both dots). The values e = ±.U mark the transitions from 
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TABLE in. Estimates for the dephasing times T2 for different values 
of temperature T and interline capacitive coupling C12. 



Dephasing time T2 (ns) 


T(K) 


C12 (aP) 




~1 10 20 50 


50 X 10"^ 


688 


588 


300 


191 


92 


150 X 10"^ 


633 


542 


275 


176 


84 


250 X 10"^ 


511 


437 


222 


142 


68 


500 X 10-^ 


306 


262 


133 


85 


41 


1 


161 


138 


70 


45 


22 


300 


0.55 


0.47 


0.24 


0.15 


0.07 



TABLE rV. Estimates of Q factors for different values of temperature 
T and inter-capacitive coupling C12. 



Q factor 





~1 10 20 50 




50 X 10^=* 


6878 


5884 


2990 


1910 


917 


150 X 10"^ 


6333 


5418 


2753 


1760 


844 


250 X 10"^ 


5108 


4369 


2220 


1418 


681 


500 X 10-=* 


3059 


2617 


1329 


850 


408 


1 


1614 


1380 


702 


448 


215 


300 


5.5 


4.7 


2.4 


1.5 


0.7 



(Ni = 1,N2 = 1) to (A^i =2,N2 = 0) or (Ni = 0, N2 
states and the breakdown of the perturbative expansion. 



2) 



The first term in Eq. ( 68 1 is quite similar to first term in Eq. 
(|6]l. Gate voltage fluctuations will couple to this spin qubit 
similarly to the case of the DQD charge qubit.^^ Therefore, 
for small biases (|£| ^ U), we can study decoherence induced 
by electromagnetic fluctuations in the spin qubit employing 
the same analysis developed in the previous sections for the 
charge qubit. We note that several studies of decoherence due 
to ot her me chanisms also present in these qubits have been 
done.niSIlIl 

The decoherence rates will depend strongly on the qubit 
operation point, given that the prefactor \dJ/de\ appearing in 
Eq. ( 68 1 varies rapidly with e. Qubit operations around this 



point require pulsing the exchange coupling J{e) for a time 
interval te, during which the qubit may be vulnerable to de- 
phasing due to electromagnetic fluctuations. An estimate of 
the corresponding decoherence time can be obtained by using 
the curve J(e) plotted in Fig. 3d of Ref. ^ Near e = —1 
mV, one finds \dJ/de\ w 10^"^. Since decoherence rates are 
proportional to the square of the bath-qubit coupling, namely, 
1 9 J/ Sep in this case, we conclude that the decoherence times 
due to electromagnetic fluctuations are 10^ larger than those 
found for the DQD charge qubits, hence ranging from tens to 
hundreds of millisecond. In practice, these times are much 
larger than te, which is typically a few hundreds of picosec- 
ond. We can conclude that gate-voltage fluctuations are also 



not a significant source of decoherence for spin qubits. 



VIII. CONCLUSIONS 

In this paper, we have modeled noise introduced by gate- 
voltage fluctuations in double-quantum-dot systems. We at- 
tempted to model the circuits leading to the DQD in a way 
that put us as close to real experimental values as possible, 
while still being able to estimate all the relevant parameters 
and calculate decoherence rates and quality factors. 

We chose to place our noise sources in our gates because 
we believe they give the largest contribution to decoherence 
during qubit operations. For additional considerations, noise 
sources could also be placed, for example, in the drain and 
source electrodes. 

We have estimated the effect of fluctuations in the elec- 
trodes feeding the quantum dots and shown the influence that 
parameters such as temperature and inter-capacitive coupling 
between electrodes have on decoherence in qubit operation. 
We have also shown that, similarly to decoherence by phonon 
coupling, temperature degrades coherence in the state super- 
positions, reinforcing the need for efficient refrigeration of the 
leads. This effect can be explained analogously to the radia- 
tion of a black body, which increases with temperature. 

Contrary to what was initially expected, it was found that a 
stronger inter-capacitive coupling between electrodes actually 
introduces stronger decoherence in the qubit system. Thus, in 
order to mitigate this effect, it is important to keep the leads 
gating each quantum dot in the system as isolated as possible 
from each other. 

There are a few possible refinements to the model presented 
here. One such improvement includes adding the electrical 
resistance in the leads, which in practice requires the use of 
a lossy transmission line model for the effective circuit. It 
may also be important to take into account the drag effect on 
the leads due to the proximity to the 2DEG. This effect will 
change the effective circuit parameters, thus influencing the 
calculation of relaxation and dephasing times. 

We have found that electromagnetic fluctuations in DQD 
systems do not introduce a dominating decoherence effect. 
The quality factors calculated for our system at room temper- 
ature (^ 210) are still well above the Q factors found in sys- 
tems under the effect of phonon coupling (~ 50).— J^— If we 
compare these results with the experimental results (~ 3 — 9) 
for Q factors, the discrepancy is even larger.l^ 

The disagreement between theoretical estimates and mea- 
sured decoherence times in charge based DQD system leads 
us to believe that there must be another possible noise source 
that accounts for the short decoherence times observed in 
these systems. For instance, it has been recently argued that 
electron-electron interactions can enhance the effect of fluctu- 
ating background charges on the charge qubits.-'^'*" 

In order to identify the leading decoherence mechanism in 
charge-based qubits, it would be very helpful if the depen- 
dence of the Q factor on the qubit operating frequency i' were 
measured. For instance, for bosonic environments, this would 
yield the spectral function. With this information in hand, one 
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could perhaps trace back the physical process underlying the 
decoherence mechanism. A candidate for such source is the 
presence of fluctuating background charges trapped in the in- 
sulating substrate or at the GaAS/GaAlAs interface. 
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